Surface-induced First Order Transition in Athermal Polymer/Nanoparticle Blends 
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We investigate the phase behavior of athermal polymer/ nanoparticle blends near a hard substrate. 
We apply the density functional theory of Tripathi and Chapman to these blends. We find a first order 
phase transition where the nanoparticles expel the polymer from the surface to form a monolayer. 
The transition density depends on the polymer length and the system bulk density. The effect is due 
to the packing entropy of the species and configurational entropy of the polymer. The simplicity of 
the system allows us to understand the so-called "entropic-push" observed in experiments. 
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Mixtures of polymers and colloidal-size particles have 
been shown to exhibit rich phase behavior, both in bulk 
and near surfaces jlj]. The corresponding behavior of 
blends of nanosize particles and polymers is still be- 
ing explored and is not well understood. The small 
size of the nanoparticles leads to both new phenomena 
J3- S 0] and potential new applications. Nanoparticles 
have been shown to migrate to interfaces such as cracks 
or substrates in experiments, and thus may allow for 
self-healing surfaces The self-assembly of nanoparti- 
cles in, for example, diblock copolymers has been stud- 
ied with the goal of producing composites with useful 
mechanical, magnetic, and optical properties |@,@]. 

However, a detailed understanding of the mecha- 
nisms by which nanoparticles interact with a polymer 
matrix near a substrate is lacking. In general, such 
systems involve a complex admixture of entropic, en- 
thalpic, and chemical interactions between the different 
components that is difficult to ferret out by experiment. 
Insight into the physics of polymer nanocomposites can 
be gained by computational studies of simplified mod- 
els in order to isolate the individual interactions. 

Additionally interesting phenomena occur even in 
relatively simple polymer nanocomposites. Recent neu- 
tron reflectivity experiments @- have found that in 
ultra-thin films of polymer /nanoparticle mixtures, the 
nanoparticles often form a monolayer on the substrate. 
Remarkably, this occurs for the case of polystyrene 
blended with polystyrene nanoparticles [8], in a regime 
where the nanoparticles are miscible in the polymer in 
the bulk |4J]. Since the polymer and nanoparticles are 
chemically compatible, the forces driving the nanoparti- 
cles to the substrate are thought to be entropic in origin 
JilfllJ. The presence of about a monolayer of nanoparti- 
cles at the substrate has the additional surprising prop- 
erty of preventing the dewetting of the films. 

In this Letter we provide theoretical evidence in sup- 
port of the hypothesis that purely entropic driving 
forces can result in nanoparticles segregating to the sub- 
strate when blended with a polymer melt. We explore 



the behavior of the most simple relevant system: a mix- 
ture of hard nanoparticles in a hard-chain polymer melt, 
near a hard surface. An understanding of this system 
is a prerequisite for understanding more complex situa- 
tions in polymer nanocomposite films (such as the effect 
of attractive interactions) since the entropic effects will 
always be present. As we will show, we find a previ- 
ously unknown surface phase transition in which the 
nanoparticles do indeed spontaneously form a mono- 
layer on the substrate. 

We employ a computationally efficient, classical den- 
sity functional theory (DFT) to explore the phase be- 
havior of our model system. DFTs can capture the 
microscopic structure and thermodynamic behaviors 
of complex fluids fill H3l . They are based on min- 
imization of a grand potential free energy functional. 
Since the minimization results in the grand potential 
free energy of the system, DFTs are well-suited to the 
study of phase behavior. Previous theoretical studies 
of polymer /nanoparticle mixtures have mostly focused 
on block copolymers and used various techniques. Re- 
cently Sides et al. used a hybrid particle-field theory 
(HPF) [13] in which the nanoparticles are treated as an 
external field and the polymer is described with self- 
consistent field theory (SCFT). Earlier work used a com- 
bined SCFT/DFT method 11411 , with the polymers de- 
scribed by SCFT and the nanoparticles by a DFT. Other 
works that treated blends with DFT include an examina- 
tion of the wetting behavior of mixtures of colloids and 
polymers near the bulk coexistence (la) , and an investi- 
gation of nanoparticles in diblock copolymer thin films 
lTl6n . Here we also treat both the nanoparticles and the 
polymer within the same, consistent theoretical frame- 
work. Our work is the first to examine particles in ho- 
mopolymer melts near a substrate with DFT. 

We use a density functional that is based on the fun- 
damental measure theory of hard sphere liquids, pio- 
neered by Rosenfeld |17||. Both the nanoparticles and the 
polymer segments are treated as hard spheres. Bonding 
constraints between the polymer segments are enforced 
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using the Wertheim-Tripathi-Chapman (WTC) bonding 
functional fl8l H^l , which is based on Wertheim's ther- 
modynamic perturbation theory (TPT1) I20I |2lll . The 
DFT is formulated in an open (/iVT) ensemble. The 
grand potential for the hard sphere /polymer blend is 

n\px(*)] = Fid\Px(*)]+hs\Pa(*)]+Fch\p*(*)] (1) 

+ Y1 dr P«( r ) [ v *(*) - M > 

where the terms on the right-hand side represent the 
Helmholtz free energies for the ideal gas, the hard 
sphere, and the chain constraints. The final term is the 
Legendre transformation where the }i a are the site chem- 
ical potentials and V a (r) is an external field. The ex- 
act form of the Rosenfeld hard sphere term, F; IS , can be 
found in | 22l] . A more accessible derivation of this func- 
tional can be found in [23]. We implemented a form of 
the chain free energy functional F c ;, (Eq. 22 from lUstl ) 
to keep track of, and solve explicitly for, segment densi- 
ties by treating each segment as a separate species. The 
Rosenfeld functional that we use is known to give very 
accurate density profiles for binary mixtures of hard 
spheres l24ll , while the WTC functional has been shown 
to accurately capture the physics of homopolymers and 
polymer blends near surfaces EES. 

Minimization of O [jOa(r)] leads to a set of nonlinear 
integral equations for the density distributions, p«(r), 
of the constituent species. We solve the DFT equa- 
tions using the Tramonto fluids DFT code j2^. Ep. E^j . 
with numerical methods as detailed elsewhere |28|]. The 
polymers consist of a chain of N freely-jointed tangent 
spheres with diameter a p while the nanoparticles are 
hard spheres with diameter a n = lap. We keep the to- 
tal packing fraction fixed at rj = 0.3665 unless other- 
wise noted, where rj = n/6{p n of 1 + PpCp), and p a and <r a 
are the species bulk density and diameter, respectively. 
The density profiles were assumed to vary in the z di- 
rection only. We performed our calculations in a large 
(L = 80(7p) box with reflective boundary conditions to 
guarantee bulk behavior in the middle of the box. The 
external field representing the hard wall is V a = 00 for 
z/ap < cr a and V a = for z/up > cr a . 

We calculated the free energy and density profiles 
of the polymer/nanoparticle blend as a function of 
nanoparticle concentration. We began by converging a 
solution at a nanoparticle density of p n Cp = 0.001. We 
then employed the arc-length continuation algorithms 
found in the LOCA software package |29l. l30|| to trace 
the free energy of the system as a function of nanoparti- 
cle density, at constant packing fraction t]. This allowed 
us to find regions of multiple solutions, with radically 
different morphologies. 

The surface free energy of a blend with N = 40 is 
shown in Fig. [TJ where the surface free energy is de- 
fined as Q s [p*(r)] = Ci\p a (r)} - Ci bulk , and Cl blllk is 



the free energy of a homogeneous bulk system with the 
same packing fraction and composition. The dark curve 
with "x" markers corresponds to the stable solutions 
of the free energy minimization for the mixture. There 
is a distinct change in slope in the curve at a density 
of jOj^Cp = 0.01263, which indicates a first order phase 
transition. The light part of the curve represents the 
metastable and unstable branches of the phase space ex- 
plored by the continuation solver. The point at which 
the curve crosses itself has two solutions with distinct 
morphologies. The dashed line represents the free en- 
ergy of a neat polymer system with the same length 
(N = 40) and packing fraction. Note that the addition 
of nanoparticles reduces the free energy of the system. 
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FIG. 1: Surface free energy versus nanoparticle density for 
N = 40 and r\ = 0.3665. The dark curve with the "x's" indi- 
cates stable configurations. The abrupt change in slope of this 
curve at p* = 0.01263 is indicative of a first order phase tran- 
sition. Solutions in the metastable and unstable regions are 
shown by the light curve. The dashed line indicates the free 
energy of a neat polymer system with N = 40 and )/ = 0.3665. 



The two coexisting density profiles found at the phase 
transition in Fig. Q] are shown in Fig. [2] A density pro- 
file converged at the coexistence density (p*) from a pro- 
file at a slightly lower nanoparticle density is shown in 
Fig. [2^. The density profile is typical for a dense liq- 
uid. We see that both the nanoparticles and polymer 
have pronounced peaks near the substrate. The poly- 
mer peak is closer to the substrate because of the smaller 
size of the polymer segments. Converging to p* t from 
a higher particle density (see Fig. Wp), we find that the 
polymer has been almost completely excluded from the 
vicinity of the substrate, and there is a large peak in 
the nanoparticle density adjacent to the substrate, indi- 
cating a large adsorption of nanoparticles. The height 
of the first peak in the nanoparticle density does not 
change significantly as we add nanoparticles above the 
transition density, so the structure of the monolayer re- 
mains the same above the transition. The contact densi- 
ties are in reasonable agreement with the pressure sum 
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rule p/kT = E a P«( z = <r*/2 + £r p /2). 
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FIG. 2: Density profiles for polymer ("x's" - left scale) and 
nanoparticles ("+'s" - right scale), with N = 40 with t] = 
0.3665 at the phase transition (p* = 0.01263). These profiles 
were computed by setting the density to p* and restarting with 
an initial guess profile which was converged at a lower (a) and 
higher (b) density. 



The excess adsorption T a is shown in Fig. [3J where T u 
is defined as T a = J °° dz (pa{z) — p«)- The vertical lines 
indicate the sudden jump in adsorption as the nanopar- 
ticle density is increased through the phase transition. 
This jump shows the polymer being expelled from the 
wall by the nanoparticles with the van der Waals loops 
being clear indicators of a first order transition. 

Thus we have found a first order phase transition in 
which the polymer is pushed away from the substrate 
and is replaced by the nanoparticles. Integrating un- 
der the first peak (0 < z/c p < 2) of the nanoparticle 
density profile in Fig. [2J) gives a density per unit area 
of 0.2866 /cjj, which corresponds to a densely packed 
monolayer on the surface with an areal coverage of 0.9. 

To our knowledge this is the first report of the cal- 
culation of an entropically-driven surface phase tran- 
sition in an athermal polymer /nanoparticle blend. We 
note that unlike other surface phase transitions in sim- 
ilar systems j3ll l32l l33Tl , our blend does not sit near a 
corresponding bulk phase transition. Previous studies 
of binary hard sphere blends have found that there is a 
fluid-fluid demixing transition only when the size ratio 
of the particles exceeds 5 : 1 l22ll . Paricaud et al. 1341 
showed (for a blend with a similar equation of state to 
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FIG. 3: Excess adsorption of blend components as a function of 
nanoparticle density. The "x's" and "+'s" indicate the polymer 
and nanoparticles, respectively. The light colored parts of the 
curves correspond to the meta/unstable branches. 



ours) that there is also a first order fluid-fluid demixing 
transition in hard polymer/particle mixtures, but only 
when the particle diameter is larger than about 5 times 
the monomer diameter. Similar immiscibility for large 
size asymmetries is predicted by Hooper and Schweizer 
l35ll . Our system should thus be in the regime of bulk 
miscibility. 

In Fig. H] we show the effects of changing the length 
of the polymer on the density of the particles p* at the 
transition. Results are shown for two systems, one with 
a solution-like packing fraction (rj = 0.3665) and the sec- 
ond with a melt-like packing fraction (rj = 0.4152). In 
both cases, p* t decreases with the length of the polymer. 
For sufficiently short chains, (N < 8 for r/ = 0.3665 and 
N < 5 for t] = 0.4152) we find no phase transition at all. 
Also, there is no phase transition for a binary mixture of 
hard spheres of diameters c p and 2c p . Thus, the phase 
transition is a polymeric effect. 

The transition must be due to the interplay of chain 
entropy and packing entropy. One driving force is the 
well-known depletion potential of a l arg e sphere im- 
mersed in a fluid of smaller spheres II36I1 . This effect 
leads to enrichment of large spheres near a substrate. 
Packing entropy should thus help to drive the nanopar- 
ticles to the substrate. We note that for particles with 
the same size as the monomer, cr„ = <7 p , a phase transi- 
tion in not found, so the particle /monomer size asym- 
metry is important. A second driving force is the loss 
of conformational entropy when the polymer is close to 
the wall. Apparently, the conformational entropy loss 
is key since the transition only happens for sufficiently 
long polymers. 

To conclude, we have identified a surface-induced 
first order phase transition in athermal poly- 
mer /nanoparticle blends. This transition is governed 
by the polymer chain length and nanoparticle con- 
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FIG. 4: Nanoparticle density at the phase transition as a func- 
tion of chain length. The "+'s" and "x's" are for t] = 0.3665 
and n = 0.4152, respectively. 



centration. The existence of a phase transition in such 
a simple system is noteworthy because the entropic 
contributions at play here will also be factors in more 
complex attractive systems. Our work takes a step to- 
wards identifying and clarifying the effects of the subtle 
interplay among entropic contributions originating 
from size anisotropy and from bonding constraints in 
polymer /nanoparticle blends. 

The existence of this surface phase transition is con- 
sistent with the observation of a monolayer of nanopar- 
ticles at the substrate in the experiments mentioned 
previously. We have demonstrated that entropic driv- 
ing forces alone are sufficient to form a monolayer 
of nanoparticles at a substrate, adding weight to the 
"entropic-push" hy pot hesis used to describe the exper- 
imental results [9J, llOll . In addition, the theory pre- 
dicts a lowered surface free energy induced by adding 
nanoparticles to a supported polymer film. The lowered 
surface free energy may contribute to the observed inhi- 
bition of dewetting caused by the addition of nanopar- 
ticles @]. 
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